


f 


j! S nN 
iby 4 ‘ ts Me , 









Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1998-06-01 


Axial heat conduction effects in laminar duct flow 


Girgin, Ibrahim 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/8020 
Copyright is reserved by the copyright owner. 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
get Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 
KNOX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 











1 $ ea 
' at fe ‘ : : sie 
yf 
Biba D844 tn prathtoty Pall eee Se 
whitey gets etd eld of ae of 
a A aS ee te ae oe Sod 
bt i Ree Geet a 





Bening 





pets tyes 
tyne aad 











ia v7 


: ee ; ' 
NPS ARC pen 
AR HIVE weed o't A art H be rare Pe 
. 1 ‘poe? pry We PR 
Li * 4 i Va F ry) ole Ou 
1998.06 LSS a Rasy 
e : ' . £ . H 
: Ane Pr ste ‘ the i I. Yer eae 
5 ‘abe * 1-4 bh Ae fis ac idlete 
GIRGIN, I ol ba Sea a 
: : yn bl. | my : Paar et : , tr * 
fg B Lee ye a NE Nage gale, wine ie alte Mest ui anise scsetne deta 
h Pay t ae state ; ey ch ani! : : yr eat of i , so tsi a pry 
\ pad h Toe payers Ly ‘ 
aera ny 
ry 






















pe 
Bt CHA UD ah gaat og #1 Pia togel 
ae ge: tat itis Sse Res berieg 
ne Wf) Hoes ee ¢) ae 4 
be Lane 18s oe wae 48 F rere 
eat Myra Mt OE Te NST POD) 
Awe dere please 

ea bade? 
2's Wa ase tas 






Pte aol 
» Sp & 































I 











2] wes ra 

mee susetieaeh eer ais, 
ty 

eS SRee ehege rite ye} Pye 


ee 

































Ome ty 


& aye ; 
8 Ween pooh at 


selne 








wut e 







































; ' 
{vate wa a tote t ped sdelgr blades 
ara sees y A aiats ay! Pay, a Piers : 
Hy} 7 we ' 2 
Hy o (st ae a4 ae hae pear) Seba ioiass ab 
Tpit A ae ab is ay Bead vhs oe eeehe eee 
oat alah oe 
re ela ooy wi) 





ese fore 
















abe ta Oa end 
- Pry * 
4 np bee ea rer Spceciads 
snes ar erae ate ara ld Wen ad 1008 5 
He ay ni Sater area 
e@gustg Gone ty ee) 


ERP ee em 


Zico 


pike Se CW 


i 
soar iytea TT a 
AStt' gs ur! a 































42'Baoht? st! 
pa 1 byt 1 
"7 ra! eh Pl ‘ ‘w 5 f PRO ECGs ae 
feta ida p 91 FP cadet yr oty te 

* Oped ping nl y les A944 anc tibet 
1 ert acid Tae vite: ‘ hh heh Ph Fait ety 

ethane fp ut “s 9 Vybsh cif vdera Ot oa say 
r) © ot b> dP aid gf 
e Bus 


9 



























mae 
Sea tgciea tant hive ii ee ° 
rvuntiipes a ety a Thee 
i a4. 
rane abated oF b Py; feed ou 
‘ Ab etsne’ 










a TOT 









uae tt 


or BF ae id a 

























> he hs rt 
as Uh ne bah ® ons: sate Fata bale ated gah rene: | 
ie te Te. Gait oo tom oA sods 





a (eis ts vee Stary 




















































































1 HS Gr gry? qt 
1p tye Pe ahs Nbr tatyled 
aay Souths? fi srefl Lak mete EAD ashy iid dat 8 see, 3 tsar A: 3 sh shadyhaparty’ 
he he s eh bpd tet Ott, ane ge’ ephaesss rvaba wily wi a amete he ey 
“ ' : 2 pre Viahate? YT Ps Fee Bs we 7} Tas ete tted fa bet et ol Aa ales wd 
4 ut t a A mail dy Vs Gals ie Situs ake + taht 4 Aw S35, APE rye eee RD ebedue, mpelne fp 
ai 4 ‘cies Ye po tase td aRe tenis! ime 0b A mt we en 
: HY ee mM} ike ¢ Sey! his aainead 1 Be 1, is be et nes, pal be We a ah aD ol oF 
; ee a gt tpat et atye oc Bahay 2 HIN te cated oe ued e vA , apie’ ta: wt ot cult Ms 2 bp site sap Pop ace! 
retae § a} ay Aq ol ei saty fatatal” CT ee .é waders 1g: 1 r DED oneal mtd? ie yh oP | grees Mabe eat line BAT 8 2 
beet Ria hided tely hatgherz aceite Renee Shs tei rie ret pr Mehra Nn TS abe techs ae 
14, aS ao ; ws r) r) own” hte e ahwlews 2 a 
regs 9 HERS” Sea Seer T eae ery nena Sp, bowed aaah Sistt Nudebatedelg oe eam: Bo Rok teas MUA ADee 
ibd (ontat Rieavisuye saad wears ayes 
*e staan eB te 11 whew al & 
Fat ake 8 rs Rgh ed and sehhsdal> 


0 fanart An hee 
4 re beh ® 4 ae 






sn aoe 
1g 









Cyto tae Pa tat eS nso hab intad os 
1 
































































































































































































































































































































































































































































































































































































































































































































































































































































































fan ¥ te ihe Say bade TERS ALY viet 
ae 3). 
a sae* atlas det otarase” ph .o ah gh wre yiirs 
nH ‘i pst eae rn) oe Aad vie tl 
: , fr 1% a E 2, 
° a ora vA : hee ¥, abog asd 2. ty 
4 Nb bq fy Ped moh? Bade? 
ide aN wet re eed Pee pap Devel el ao 
AS; ite AOE» ; fas Parade a Sheep ¥ . c ‘arab. Lt ee AY 
* ° ty ‘ae \ 5] faite’ a8, gee ble Boa De' . 
. adgtider ts a eas . Bx i shit Pa eye ae se anne acdcthe nhg spt abets rite shel Amid gate het 
% ObRbs abe an ouas ARE abel er rp ter rabeh tee F.8 iy oe Nh st RAS Gated SoS 04 0 mba ebatve til a) Sa. 6h. 
‘ 1 4 b " pyeard ts Aah Sabhok METI OF) Pe) eg ive Baieatl (F:Ny Bok. gheha® Sede tenants at wDeecee as 96 te rod Se 
iets too. il : F "1 iG rq bia thep Raby Soy Wa ms i” } ey) © nf oviy ty Lotasd hb ut oe Rsceepity bo} iw bale Asal 8 pao se bet e%s 
; AL , eg eh 3p " 2i4at ye a a 133 st na olen fe salen’ jn He Owe ie satare hb jo Maw Pi ! 
7 er iain ams | tee spo age ate Pte aa ci i Recaaiin osaee a 
eon | rite Share 4 i 
1. an ; i » Pad aan gee 
1 ) st - ‘ ~ i ee 2168 
' mt oy ¢ » 
‘ - J tv ' : R 
1 ; 1 P 1: ‘ 4 ' ¢ ; Ny eaaiet aig > Th “idyrete Zag! ee fy ene be pease Fo ds 
‘ ae ‘ ' ‘ Ae ' - ad ad 46. -“/? 
aaa cn Bd, fogizet a 8 <i ign 18 AN, ni adage pelts! ASe “at digecoae 
wD ~ oa ' s 3 not ' $4 ‘ ; ooh ARG - 
ee Wr ' ‘ ’ es 2 . ond su asia: oe ~ Be 
4 ee ba ups 2 
3 aa we el ' " H ' Bae ay legge! 
« 1 ? ' ' . : 
1 1 » af 1 - vk a & 
« u& ¥ ‘ { 
: 5 pai ’ ATT Ae may 
san ' 44 4 saan si4 Nas 4B. Ba AUS Wee wey ChEMS : 
! 1 ® > PPhed es Withers 9) 4 a" = [oe Ante 
no : bute i £52 a en 
" +6 1 . ' ’ -* 7) debe a wae" 
; 6 ate 5 te. SS iieaces at thm of 
, eas ; We hate a 
a : peas Pre wnt ab 809 860 Sab bLietated : ane Tat wwe” how” 28 
5 : 1 “28 athe ui : rc atind dnt Be laAstahe tas 34, mated 2 eee ates 
: ; weit oN + rest ‘ eet Kuen ir PTA Poe Le Pate? ys 1 tae em ean Nee Seo Nl 
S ve ¢" : ; : Ta ee €igey ; He tetab: al eve ‘tat aha’ Rasa tate teh ii cws? Sw te Ae 
f ' aC a ; ry : Fs rang Vela ; : tte inesteieae ahs 74 . é 
viet ; 1 Magee gem hs ‘ ' Ser oi See aaa 
J 1 4 mae ut "4 t ‘ " eS t 
: ae Tee : 
, ' hy 5 WE Fr tary 
: 4 ‘yt r et Ss gin scts 
we ' ‘ ae bee? , 
» 7) t * 1 % : 
1 ‘ o's .) a4 
‘ ' t t 4 : urd 
4 , ia a Bs se oh 
« ‘ ae a. as be aS ate 
F ’ ' mt ' 4 3 Mr € 
’ . 1¢f 
Piaras sone vod ame neve 
’ ’ ne a st ee 
v f Oa f 
‘ J ry 
bd r 
\ ' ; ‘ : eee - 
ba a ‘i “, 
, ? . Z t & f d : chet 4 ise ope Bette ayers Ae 
4 ® 4 ’ of Metaiwighy ott & 
" ) ; ‘ ae H ns * Wp ‘ he oe 
i F re 
, ‘ oy HY ’ 
te , we + ( 
teat . 1 if a ‘. 8 4 ine A 
t 2 4 ‘ . ‘ eee if gts Raerte ta) 
‘ : ta sh Fastrt ss 
: Pag ti hie. | © pt 1 
{es s ae 1 P t ¢ 
ee ore 
! é ‘ ; 
A ; i, 4 
‘ a aig a 2 : 5 meh 
‘ F - a 4 he | 
’ 1 $¢° ca. ne v 
’ r7 sa ‘ 
2 ’ % ’ q 
| 9 
¢ @ é a t £ te wu 3 
Soe 3 e. ae? ae nm | a 
3 ; 
‘ nm ve eweu ‘ . reenter” 
‘ , | : 
‘ : so. macte re ooo 8 
: . pun ye o id . ores 
Lue SE ce vant weezy te 
; i > oF ‘ 4 a Y x i" ; ae te 
sae ih dg" rhue uy 
fe 
: ‘ 
q os ® F ri See 4 
4 o 
F 770 
‘ + yw trys, % 
J Z B 
" 1 H = HE é 
# : vir peat ly Sar nyeene “pets mg ers 
ae iene tartare Scamto riick 
ata i bd ¥r fg > rf Ly sa 
i , 1 v4 rere i BEB ree hr 4 Age: es 
ar Uo age AOE ‘ 
Bs $e : © nf 
sous : ge 
4 Moen dtd te“? se } Pe a tel ie Piotr oo ad 
a , ' arc Lae th as “ef by nan Ant ito tevt tds qeigte? es a gor8 (tart Merer é PPT) 
. , she tadatalgt aes ype steel er ere il el ule op Uals Sh wou a eits Sere re wens 
; i He ye ess 0A Alle is a8 gt qoere heer Pes das fasgeore Prt igs yt Teed" " ae BA Wie ; “ote we! 
io Se alae are ae Bae tree wheter es Taka eH OT na ‘ ree ots SPP OPS at of el 2) te eke date ae ¢ 
al .d " : coeit| ee hakange [Ue ee J Fe ns tay, ee or tee eas yeas rie? Bae Finer iat Teh eX 
é : Agh ge GU Et 31 ¥ oe orepletig te A Hs (SiS H EwHiee Fete Ltad eed ts! i886 gagtere tag 
SAS Oi ee eetaatives segures* yeas rarer’ a 7S steer erin inssea* gy (gh (aN 20 «*e srr gusenat: 
ny ee eee Lee ey, ats pig cpseng tree Se trae RE gE 49% useage fi areP eee ey dot Poe etn iy 
ane atte $4 sreeays uF reer aen Chey VE sa ys rd ol be Oe A orcee the Goes Ae FOL ee WG SY STS 
aa gets 105 of | Zot QOTUMHTN Faes ; pb pastes AF pra het ge Beet ra at oat mg dames a8 ya: 
mune L : pepe een i pd be? 4 pe? ates erer re eat ae = att de ¢ nawrgy ae 
aft : ant the dor doers? meee LA Beit 2.4; i Frau Hi ig 44 otha Pp cata tn tae er ae & rot? i mete tin 
re a * og peptasett y ser ic atytigt goa satetyandet? b cfr : PW aretha blader ie eee HE ie eye gle iat st 
SCC) OL at Ay SO ie PI ie otto bts Be pe bape A me POP Dy var 7 Teeter eter eete Dh ekes OE MI SN eee Dastteresne 
mr ie ‘ot oeemee 8 F Ea is TU ES ab sd hac odes wt 18k "92 peed gt gteerhh ro prot Hae ot 89st ed Fo 4 340g tg ga i* Peerter ae POT TP Tag 
cok ele : oe eres 2, ogee ‘ ius ee fragt ge IR EEA eS oer firs RN Rt ee te Aa 4 tg? CONSETT oak ile -angrayen oe gti Mgt ess! 
2 AE ae et)! aren : ae H LTS ate Pad eeaty® poe re abten By alt of aed pryhietes geet pals tigt yar ree ron ed ote, ee He PT DAF be dll oe satgnens“ceeee ess 
Py ae 2 CM Bie re Hae y Sagbedg ret 490% “aagr th. § we sree’ ep e THM Me! Seaatal Phe eared 20a" stelh etd 7 er ay bake tm ete! De -B-| per 
pr yegeee yo Be He rT TL dap an 1< P| eset aftss = Ft yaPhe th Up Age pent ar ea (age af usp 99 Ue qeevisl e 7 PRE ey aa te all eles et yt OF pac gt; Simic rcee ae stage aed 
: cin ee ae Giese yas CT ad 4 ee heh rand of Li La bt al af ated Dotti nF eee Rms Page sinre ya pt wer: ae 40 145 fet phe ga eka. Ae LA 8 gee tg Bree? eeceles Cu? Rt atid 
: Mg tO os oogejeu7l te gas ¥ NG pet ota pate gear C14 SSivierer riparee Pea titel 4a" HIM ie cen pat eet a eager eaten rust ate 
* rar wot i Call Micaela aa oh lay th niy® ee Lert? gases hy eure! Patat ieee ‘ He rs sedotys jepeeyert Hol gl of ferret hb bd Sf a4 of 
: ae ee oC ae ne emypse Atte eb ides paged Foi ea tLe diye Beneckans a atcha leh aoe eset mate tongs Let f apse? : tates rer 
‘a eet Ry asgeye ® Lead yeotd ome limes ’ uty : 5 iyhete Sind gn hetee ge Igy roetet Fae Sp panes Petar chy 24 Jee Pa we SF ah ne 
ry a4 | vat coete attest? ‘ an ¢@ grate ot itp ais u digs rseeae ct Hes” £ Gh t Part! 72" ser a aia testards errs ye el ak 2 pe EY ie tart, 8 
' ry atte oe er a Fa = 10 any ag! 2a teteer) +8 same ett age au te" rhe <i) reese op vr yeah pg ope by a . es [ee et Ap re 
ve t af a > 4 ptr err ajh pert era* 494 08 ve 5 a 
ae 1 a i paler meer aeaver Gepaes 1 cs fy f s we?  # Aas S| | Lrahytghe nif appt seat eter sease'y oe a A 7 pei hice afr page sed? foc OP ergs ie 
‘| ra 1 gp gies Satan gp agere eng yates. Af "yhceetete but Oe lips ah pe ah A OS sande amy! oye a Cie 0 ttalet Gr tr egta@ WANTS ite qe taroniane rete. Serer ares nen ae TS - 
40 vem ‘ eee prema s pa(sae tase owed Pry Mere ir en Ciecererag fom a? geal * UTA a ip ads FE A esta e) ee ” Sak el ch 4 Se ae Paes whois Ca thls dete 
i] i4e Ate Hs ese pornne cedar lor a rans thn ALA a & Sop 5 at oe ge fe eee pet Ie eeh ye PACES 
per ie rat) faghent ehpuisinsagy Oars + pipe RTC det hed bk eh 6 Tagryeses 0a? Z es el een HO Ly pee Ong gmy rte ee ty 
meet Mee ee iyecetagace gee ieee F431 SoG haseay esate Te Tt Seliicmtgsrenceusinat pre Clonee A Ara roms etr lL (See | OEE NS 
rh OS dei sare bye coe pee rey tates pees gty bebe ie Mri Pn Peace eter abe ed rn rir ooeurse Meine MEN TALL, pete Has 1s abe be $Y 
tees 4 ade as phaeed jel fyagassite'y ‘ tedaeus cphatgae yb i Page rae wahitss ere er: 1 Pn eee Corps AP ETN OLIN Se 
J . Sone Lapel fe eae 4 » ; Bg @oye oyoyea dear ey cata’ +44 agent a SEAT deeded eT ate Mant oa Ang “ pe) 
poe = ouvedacyeeta gas 38 y oe sey 48 fT ‘aie gids ja ca eg aoc nses 43 vie eer ign styo8 seta ceal cto rites. sea ne nraee beter getinet rte ; 
auth by gtycatles ent Yen es eigen rae ees 285 es ete ey fetyd sen fue rrr pel did pt Bagh pee Sosqnetarn: Hapeyegver crete tars? raneg Y 
‘ en at esas Ore e He ron] > nts pate Me ne rir ¥“¢ elt eh aawer cat Manet nt ey hig? OF” eA wee Ageo"e epee pea fq uO 58 94 ne ti eee 
Chari iT a Cay POE CO Utes Bid | jee Het ee pees ab eT OI gig eas Fares wer oP Peer ozpree” wiseacaty Lor teh: et er ree th Anat 
ater) ay ire 4a tisthogr 4g! tr ore.d's had BT ~ eras 440° aetit ap Sher “betreess FLL pthe bes toh bode phe z, 
ar eee fe 1) er Arsene sf aia? et" mine 14% teres Ls areety: yegtses® Lr seein cA ew oe a tate Ftd Pers Cees 
afar a Eat ge y pe peeith ee! ate a ats 18 we eietay at : anges erg wayre taste wrgiree mnie pases en é red FT et ale 9 
Th of ty a a a 39nd ae Heya tse & Asam ere ris ees pe ry padre CP ati) a Srshcdoles F etetg® ate ri aqee? om * a ot oy Ph eg ma oa 
al ol LP tear PMN RY ay 9 ba era rtd Pl vjure {Hoey jeeen ee hays W ogee onety re rp eee bos Sm Mp re 
err ig) CODE TRL Le oe eiyts wat Mies ' i Nand ee Seen is BT EAS ne on pe cmene NIN 
fi : ‘plbncell ere of cea led aggre Pew tgtiapesasgey! Se eC OL OOS ON oper rine 
rete Pom Ai fe rea at em ranttne Scant ect ¢ a pal ae foe nee Ore r®, 
MJ el) athe hehe i F Ota ieshed es maa tft ie Lappe 4, : ne aa bef $ @ elesiy egey geo f ” ae" no é a ? } 
2 - a al ah 4 1 TD dal dad, 4b ene ey oh eel a 1¥ ¥ ePAypeerat cary ® tt ed Ae Pat ke, eget opne ep gee a 
eae eet! teat me iy naa aieeee FE dawail Pay Ae rhea made’ sages ppeteer er maeant ie ere sree or seo vr gbaeronartesee'e. a bd thera eater emer wee 
do oatete att etl ad stark poe ¢ tog Mata MVEA ELL® OA egtae pe eteweradad ( ecilt 1° 284 8ue ee e arn bs Road : id “re ‘stp ld eet 
a8 ’ he apy orate taphe 9 1tyes dn® t4gtatw wave ras epg a eh. 0d 90 ng my rate FITTS ey a an So etl 
septal Pise : IN on 18 iho Sept rtee bite phy Haye piel enre qeicveuest, : fa acaba rate ines sy" Lat 92 : z: (ee ad oe a8 ca oon Rc eae) 
aon ys Tent wb eens PD ee Pid ESotaea we i484 to damveree DE Enpreeh ai? of erin rag ak BL ae nd tal reasasee e Tene Te : 5 at ae ry eae e nese <7 ° et Gf fr 
Is on Ap deeny ’ *% pt tetvn peal be as urge ered yah Sapo orga nat* Beas ergs rerigace ey eye AM eee “ peat rege Siti ces 
; tt tipioa wyrat esate ay ard d Bache rerpstchtensw tt meth teeta La vi 1 tae forces phase Fs 'F [Aebietee Lp Patong 4 pts qe RT 
ah aes? ri) I eee Pre Hace ake tae) eases 14 Gees aye Hone ’ Pier tet eet cman e} weaeeeahs prs wa sengt aN! fete ee 
f = ig t : revel ante Noted magrle tao wages ad Ae ay 3 sites het fake seotig: ay ont vans RAR ert pial ar 
“ > so94 oi e@ . tee 5 * 4 A} oye: gare s* ” oo “ 
i ol bitte CaM itt aie be Poni ph ef beat ea ft ae + een aalaoe tee reba Sa Sias ei Ohes Sits : creer 
uh Pans baie ve ign yh gelled dean Se 02a Ke ots _~ nies ” 
ne Mees “ee ee Ses Ui ah sti tet dase ite eae one a ara Sages eeruret 
rij a achat. : wh? 849s he yn tees ae 10 Septet + es re A Tag ee at: eee ohare coeenere ra aogr ret 
' ‘a euyouees e2af epee} 4,9 a a " e : etertet " peu tg tes O8e" A ea het bathk eke dhe oe : 
besa waghees ines year: Be pe ode pays ape eh NEES AY Te ge4 F Peart ye tial tad sotnig is iatator Py “ weiter pad nes tpt e La tees Rema eT shee eck oe weed me Rte ay eee 
tinerree Be dinar pir aw: oe sleeps pal eee at Pe, ae cae mee pgatys mde gs! oe sateen. RT ST ha Basel Le yo sped a Some? ISA er ECT 
iP Leeneves Cpe dae s iret hia 4t08 ater ner eae SGT eb onea a caaacs feu catig Maly 0614 fa be oh oie ie a gadabh oe heaped bet bf Ne cals eh smth th ik Lari) 
x co bats Sige Ad eee ey aeee ‘pusevst ante 2a? ty event P4008 rats anon bie Podgt atts eateageatatee ee ee ig ast morgane OIE Fee y wee PUNT WE 
a hel ear bates ° iy A aut eure’ gonad ep oy e agg Ws aynteed roetaTetee te oriy=e? Ieee aS 
HA ectnte Cal ciad kerf tenn ML oa ipl diane cater, rhea minigun a tga ecg oe eater 
‘ rat 4 os 4 eo a 
roe ase ONE fe iotihet cM esc aet oe Pama » : 4 a gee estenates cates : bade sont Re ete naancyare=s-owne€ 
rt Fe Li ad Tee " I! Ie pahe F1Ws 5 eat pto-ae sine pon oats OE De Me™ pw western lt terre PP ek tty Pu bees Leet oe oie hand OC WTR ER 2 Dy Pate OH ® “ el ae pt led ia-<f ere rq a aeee vootdtee ne oe rte” 
aa eennn vt ate fe Mee crcel ead Pe! pee age evi an yauget ewee | fy wth 5 9 07 ies! platen somone umebarit a pereite Ley? arene! a bd a 8 ye eno Ne tbtetel 9! 
, GO v98) ha aeepep se ee ea spat veces tas tga cane o¢ [ost do 47 8s | oes rt 4 ¢ 5 ’ bape rs fo Poh -c-ro-he E 
yee 1saetes 014 r a= as pibe , Jest rede aeeth 8 seesd ed CHW h es ue Her eT td tas aie ten ipaup betas: grat gear mae 26 oh 4 1 8Ge< Ow EPO! one ee tk pd 
eat, Libs bee a slieamerten ee ea het pee DAs Fee Neen TYE mie Parente Latest ey ie eats ware nlerads ; ha oh ae i lieth aan Penge Tae TOES POO: 
Gee a te ff ghetes i ’ series {ti A ate UP hee itd ‘a al oneneute? pale! cepts ab he 4 ” wut ZANOTTI R ME AES Ue eS te a ROUke. “9 0 orme. ee ti-88 Pe eas ide 
ied) oe cee tg! Hf 4 ne Pe caneett fe fees dy uticiigenree G4 imstte he iho 0) Bas themnctne set arbor cae tae om Petter ee 
* Reale tae Sas yes ¢ Th a ae x bre’ tbe ram sore et a 
otf ee  aenese 1 pies sees itoader ae ’ i actemiericaag Haines PLA a mincing tagts 2 oa" ce eeeae eaten meereee eS 
ease C let tay ony eit shes " yes vate eg PR ipererrs ye Tepe arieaerecc™ orarenesine ria tasetd oetone > Og@D' 
lee NE de ate a tase Fitedot weiss / pie ite . Seat te fe atagy : creep eet Fo Pars y Ere age (tee ree ae hat awe rte 
fe ned Yee Spe 4s qi ciewbeier arse be 70ke 1 oF iW = aie elttuntg ijeal t ing resort °8 sare * ATOR Ai as a a Oe OT eal Spain poe sme 
yee te el Mate ee de pos Hee! ST able anion: r) a goes baie adets ay & Th ages me pproet gen gies cree? Sys lapta ceey rose 
Ll eras 4 Wy fen eae ee 18 AD eseph C1 feed 1 aso Oley te ” maaegey tte, eg es , ae : wa, vent es edt oh 7 tem rh 4 eddies Terese ce come? ceytarg sence ooesst’ re 
a ie oabateati pansy 4h outs 7 Sh Eee Pes fe naa sup ies at P ertserrtw and Sah Tafel : ia ator wart ste Sac p A eth Satta tata cen yawns 
See pete # Tht tae greed ah Bg 08 i at gee ads et veute sea fecatene att "v Mier eee df dey ie oe es tatees Boe i i tinea 2, ‘ St ke ante beet riiataern td 8‘ are Por pov ate" searoce rose ag ey Freee E ang 
ee area aH ia Port il ri feud he 30 yeocedsense bad oe ent te gisyrcns y rrasey “unesecduesyedye a! qierrest oe epee cy tat parcmeete perteaten sgt Key hae - 
| mis itetett bet ete wet ratak rity rasan cane Fgact eauty firme lg 492” 2? job th FLD a retrlreetstat ep 699 9 2g oH M9’, 9.39 68 seieares sas Pi ald ep a tpg or 
tna <del ed 38 or avid ra yo} otis wae fy. 1e'? "xe ree ne ys vine Jape sé yteale, ‘fi aecciepi get ss ent ia nett a beaced - er 
fot Ay b gel ‘ ; aaa ofA > atl eter sige? ‘ Fi 
Ls eye $ arog? 4 res + pe fe nig hee 1 Wi a ced nt ie 20-9 ye Te 1 Ste rattan peices 
pains £2 plat Nh. ae paeeinta res wel’ ‘Ie 4 meas ty ¢ 
saeaeree be ete oxy 20 eae sp Pears? GEE NMC, igzt* b + 3 ‘ \ ry oa ey i ee sand earsraresd 
Pip sencuenetieh eis pierasonty S14 81d past He Bape Mes; nn ain aerate at oh Satie pe ale 4 eye 
seridyaesent ie aes Latha nati pole PTT PT TT ahh ied rp at Hate renee rf isl ped ra ta #0. + Har 
: or Z 94 weed esha? ‘ is eee et aticdan te ofr: vibe aeaehie ae ‘4 ee i “a4. Te ihe ari navel ssenter 
te A Ae fadiunee PRON Peale ett STR te ns pe bi if y Be wap ateaie We ses grt eyed tiene 
ior eategectewe eet aa av at “ey ait aye ote bk pate U ; ous t ase ach pt 
Su} ane 1 ne Safi medert tf is nana ip isasetts Gy i rt y 
; patie MU laeauns ates feof vhs Sogo lta eatucei ad 
wets ow opt hae? leet f , 
PHM ou da ape tun Es i put ae a4. ty Pees 
" eh tet 4 tod poted ale iy Ad Pe as = 
tas Ahi 
1144 
tad 


'» oigrete ec ca Ses 
A nadig 
tr 


os a 






’ foriwe 
iy Past ARAL #0 






tuts gat oy 
is a NOt es id, pate rls hy 
1 sel oe orepuly fate . 
ma Mere fe pak eh pened pea ob at ees dadee anes fu te OT oul 
a y het hid ona tt 
hs i eria 4 7) 
oot 


DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5101 








NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





THESIS 


AXIAL CONDUCTION EFFECTS IN LAMINAR 
DUCT FLOWS 


by 


Ibrahim Girgin 


June, 1998 


Thesis Advisor: Ashok Gopinath 





Approved for public release; distribution is unlimited. 














REPORT DOCUMENTATION PAGE Fora Sperce ccs ONG Ne 000i 186 


Public reporung burden for this collecuon of information is esumated to average | hour per response, including the ume for reviewing instrucuon, searching existing data sources, 

gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 

collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 

| Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) Washington DC 
20503. 
















i AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
June 1998 Master’s Thesis 
4, TITLE AND SUBTITLE a: FUNDING NUMBERS 
AXIAL CONDUCTION EFFECTS IN LAMINAR DUCT FLOWS 


6. AUTHOR(S) Girgin,Ibrahim 





7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) . PERFORMING 
Naval Postgraduate School Sesh. ea 
Monterey, CA 93943-5000 PORT NUMBER 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
Research and Development Office AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 


12a. DISTRIBUTION/AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 
Approved for public release; distribution is unlimited. 


13. ABSTRACT (maximum 200 words) 

A numerical model for heat transfer in laminar duct flows has been developed using the finite difference method to 
explore the significance and extent of “back-conduction” at low Peclet numbers. The calculations have been carried out for 
flows between parallel plates and in circular tubes by using different Peclet numbers in the range of 0.05 to 100. For both 
situations constant heat flux and constant wall temperature boundary conditions were used. The validity of the results has 
been checked by comparison with some existing results in the literature, and extended to a wider range of parameters 
including conjugate wall conduction effects. The results are presented for bulk mean temperature variation, Nusselt number 
behavior, and energy absorbed before the heated section, for cases with and without wall conduction. Such axial conduction 
effects may be an important feature in the thermal characterization of microtubes, which are to be used in microheat 
exchangers. 





14. SUBJECT TERMS 





16. PRICE CODE 


17. SECURITY CLASSIFICA- | 18. SECURITY CLASSIFI- 19. SECURITY CLASSIFICA- | 20. LIMITATION OF 
TION OF REPORT CATION OF THIS PAGE TION OF ABSTRACT ABSTRACT 


Unclassified Unclassified Unclassified UL 





NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 298-102 





DUDLEY KNO | 

X [ lB A ne 
NAVAL POSTGRapi > 
MONTERE 


Approved for public release; distribution is unlimited. 


AXIAL HEAT CONDUCTION EFFECTS IN 
LAMINAR DUCT FLOWS 


Ibrahim Girgin 
Lieutenant Junior Grade, Turkish Navy 
B.S., Turkish Naval Academy, 1992 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
June 1998 


UATE Sc 


YOL 





ABSTRACT 


A numerical model for heat transfer in laminar duct flows has been developed 
using the finite difference method to explore the significance and extent of “back- 
conduction” at low Peclet numbers. The calculations have been carried out for flows 
between parallel plates and in circular tubes by using different Peclet numbers in the 
range of 0.05 to 100. For both situations constant heat flux and constant wall temperature 
boundary conditions were used. The validity of the results has been checked by 
comparison with some existing results in the literature, and extend to a wider range of 
parameters including conjugate wall conduction effects. The results are presented for bulk 
mean temperature variation, Nusselt number behavior, and energy absorbed before the 
heated section, for cases with and without wall conduction. Such axial conduction effects 
may be an important feature in the thermal characterization of microtubes, which are to 


be used in microheat exchangers. 
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I. INTRODUCTION 


The heat transfer problem of laminar fluid flow in ducts which is known as the 
Graetz problem has many applications in technology and has been studied extensively 
since Graetz (1885). The classical Graetz problem considers the forced convection heat 
transfer of the fluids flowing in ducts neglecting axial conduction effects in the fluid or 
the wall. 

This study is on the axial conduction effects of the flows whose Peclet numbers 
range is from those of liquid metals where axial conduction may not be neglected, to 
those of oils where the axial conduction has almost no effect on the temperature 
distribution inside the ducts. Developing micromachines nowadays and in the future, the 
axial conduction effects may be a very important feature in microheat exchangers. 

The effects of axial conduction on the Hagen-Poiseuille flows between parallel plates and 
in circular cylinders for constant wall temperature and constant heat flux have been 
studied in this thesis. The governing equations have been solved using the finite 
difference numerical scheme. The validity of the results has been checked by comparison 
with some existing results in the literature and extended to a wider range of parameters 
including conjugate wall conduction effects. The results are presented for bulk mean 
temperature variation, Nusselt number behavior, and energy absorbed before the heated 


section, for cases with and without wall conduction. 





I. BACKGROUND 


A. HEAT TRANSFER IN LAMINAR DUCT FLOWS 


Consider a flow in a circular tube where fluid enters the tube with a uniform 
velocity. The viscous effects are important in this flow and a boundary layer develops as 
x increases. The flow is “fully developed” in the region where the velocity gradient 


doesn’t change anymore with increasing x (Figure 1). 


Pi Heat transfer by way of constant heat 
flux or constant wall temperature 





Figure 1. Laminar fully developed flow 


The Reynolds Number for the duct flows is defined as 
VD 
Rey =— (oa) 
Vv 


where V is the mean fluid velocity, v is kinematic viscosity, v = a and D,, is the 
fo) 


hydraulic diameter, p _ 4,Cress_Sectional_Area, For the laminar duct flows, 
q Wetted Perimeter 


Rep < 2300. 
The Prandtl] number relates the temperature distribution to the velocity distribution 


and defined as 


(2.2) 


where ais the thermal diffusivity . If the Prandtl number is one, the velocity and the 
temperature profiles develop together and at the same rate. The Prandtl number is: 

Pr<<1___ for the liquid metals 

Pr= 1 __ forthe gases 

Pr>>1 for the oils 

For the liquid metals, where the Pr<<1, the energy diffusion rate is much more 
than the momentum diffusion rate. It is opposite for the oils, in which Pr>>1 and the 
velocity profile develops faster than the temperature profile in this case. The value of the 
Prandtl number strongly affects the relative growth of the velocity and thermal boundary 
layers. 

Peclet number is defined as 

Pe-Keala (2.3) 

and it is a measure of the quantity of the axial heat conduction effects in the fluid. The 
axial conduction is assumed negligible for Pe>10 and the axial conduction term can be 
assumed small in the governing equation for this case. But the axial conduction effects 
can be significant when the Peclet number is smaller. The purpose of this numerical study 
is to show the importance of the axial conduction effects and the heat absorption upto the 
entrance from where the heating starts when the Peclet number is small. 


The Nusselt number is defined as 

(2.4) 
where k is the thermal conductivity of the fluid. It is a non-dimensional number where the 
convection coefficient h is calculated. 


The laminar duct flow being considered is assumed to be hydrodynamically fully 


developed before any heating effects are considered. The velocity profile in the laminar 


u(r) = avia-(=] (2.5) 


fully developed region 1s: 


For the circular tube: 


For the parallel plates: 


u(y) = fy) — (2) (2.6) 


where r and y are measured from the centerline, ry is tube radius and h is the distance 
from the centerline to the plates. [Ref.1] 

1. The Heat Transfer 

The heat transfer for a duct flow can be expressed using Newton’s law of cooling, 
q =h(Tw-Tm), where h is the convection heat transfer coefficient, T,, is the wall 


temperature and Ty, is the mean temperature of the fluid, where Ty is 


Jouc,Tda 
T, =e (20) 


mc, 
For the constant cy and incompressible flow through the circular tube, T,, is 


y To 
T, => [utrdr (2.8) 
0 


Oo 





where for the flow between the parallel plates 1s 


m 


i 
T =— Akos 2.9 
al y (2.9) 


y and r are measured from the centerline, r, is the radius of the tube and h is the distance 


from the centerline to the plate. [Ref. 1] 


The energy equations for the cylindrical and rectangular coordinates are 


fy es aura eer pee +2(e2) = 0 (2.10) 
Ox OF) FOF Oh) “aC Ol mOOm @ OG Ox 


Oi Oi Ci | O(, Ot Cae or Col 
DY 0 eet AY ee Sf Oe Ve || | he ee ee |) 
Ox Oy OZ NOX. "OX CV Cy) Oz “O27 


and 


? 


where “7” is the enthalpy. These equations are applied to the problem including the 
boundary conditions and the temperature distribution is calculated analytically or 
numerically. [Ref. 2] 

2. The Thermally Fully Developed Flow 

In thermally fully developed case, the relative shape of the profile no longer 
changes although the temperature profile changes with x. The criterion for a fully 
developed temperature profile for a circular tube is 


ral sana | =const (2.12) 
Or\ TI -T =, 


m 





This condition is reached in a duct flow whether the case is uniform heat flux (q_ 
through the wall is constant) or uniform wall temperature (The wall temperature is 
constant). These two cases have a lot of applications in engineering. For example, 
constant wall temperature in boiling or condensation, or constant heat flux using the 
electrical heater. In the thermally fully developed flow, the convection coefficient “h” is 


constant, independent of x. [Ref. 1] 


B. BACKGROUND STUDIES 

The heat transfer solution for laminar fully developed parabolic velocity profile 
flow inside a circular tube in the thermal region and subject to uniform tube wall 
temperature was treated for the first time by Graetz in 1883 and is known in the heat 
transfer literature as Graetz problem. In his problem Graetz neglected the axial 
conduction and solved the problem. [Ref. 3] 


The problem statement is: 








2 
a(-2)S el (2.13) 
2 One Ora te OF, 


8, =0 ate, =1 


aa, ate = 0 
Ge 
G4— | atte | 
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where 
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DoS = Re, Pr 

and To is the wall temperature, Tj, is the initial temperature, r is the radial distance, x is 
the axial distance, D is the diameter, Rep is the Reynolds number and Pr is the Prandtl 
number. [Ref. 4] 

The velocity profile is fully developed much before the temperature profile when 
the Prandtl number is high relative to one. For such situations, Graetz solution is 
reasonably well justified. But the Prandtl number is very small when the fluid is liquid 
metal and therefore the axial conduction effects may not be neglected. For such cases, 
since the axial conduction effects are not included to the result, the solution may give 
inaccurate results. 

Michelsen and Villadsen investigated the Graetz problem with axial heat 
conduction for circular tube constant wall temperature case by using a numerical 
procedure. They used a method that is the combination of orthogonal collocation and 
matrix diagonalization. They didn’t include the wall conduction effects in the problem. 


The problem may be defined as: 














0-0) 21S) = ts (2.14) 
Oy xXOx\ Ox) Pe’ Oy 
where 
2R(V . =a 
Pe= (Ve)PCp os eae = — Jy ie Ty 
k R PeR T= gue 


The boundary conditions are: 


yoo d=) 


Ox 

Le tor ea) 
Xa Ox 

0=0 fOr = 20 


where V, is the average fluid velocity, R is the tube radius, z is the axial distance, r is the 


radial distance, Ty is the fluid temperature at z > -oo, T is the bulk mean temperature 


and To is the wall temperature at z > 0. They plotted the results, Nusselt number vs. y, 


heat flux vs. y and @ vs. y where the axial heat conduction effects are easily seen. But 
they didn’t include the axial heat transfer inside the wall, which increases the fluid 
temperature very much upto the entrance as the wall conductivity gets larger. [Ref. 5] 

X. Yin and H. H. Bau included the axial conduction effects of the wall to the 
internal flow through circular tube, by using two parameters, duct’s outer/inner radius 
ratio and fluid/wall thermal conductivity ratio. They plotted the graphs of “temperature 
vs. radius of the circular tube” and “Nusselt number vs. Peclet number’, including the 


axial fluid and wall conduction effects. [Ref. 6] 


C. NUMERICAL METHOD 

The governing equations are steady state elliptic partial differential equations. The 
temperature T(x,y) throughout the domain must satisfy the governing equation and the 
boundary conditions along the entire boundary. The finite difference method was used in 
this numerical study to express the governing equations numerically. The central, forward 
and backward difference expressions used in this numerical study for the grid points in 


figure(2) are as follows : 


LCi, j+1) 


li-j) ayy ity) 


‘i, j=) 





Figure (2). The grid points in finite difference methods 


The central-difference expressions for the first and second derivatives with error 


order of h’ are: 


Ot _ Misty) 7 -) 
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The forward-difference expression with error order of h’ is: 
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: : : oie 
The backward-difference expression with error order of h* 1s: 


Ot on Ei, +2) =F At, f+) = 3b, 3) [Ref. 7] 


ay 2h, 


(2.15a) 


(2.15b) 


(2.16a) 


(2.16b) 


(2.17) 


(2.18) 


Substituting these finite-difference expressions into the governing equations, the problem 


is reduced to a set of linear algebraic equations, depending on the number of grid points 


used. Usually large number of grid points is desirable, but the number of equations to be 


solved becomes too large in this case. Gauss Seidel iteration method can be used to obtain 


9 


a solution. All the temperatures of the grid points are then make equal to the initial 
temperature at the beginning of the numerical solution. The temperatures at the grid 
points are calculated by using the governing equations and the boundary conditions. After 
the temperature of a grid point is calculated using the finite difference equations, this 
updated temperature is used for the next grid and the temperature values converge to the 
exact temperature a little more in every iteration. A technique called over relaxation is 
used to speed the convergence of the Gauss Seidel method when applied to the elliptic 


partial differential equations. This technique 1s: 


ty new) = Way new +U—Wia ys ola) (242) 


w is known as the relaxation parameter. For overrelaxation w is between 1 and 2; for 
underrelaxation between 0 and 1. For this method, overrelaxation must be used to speed 
up the convergence. [Ref. 7] 

Forsythe and Wasow show that for a 45 by 45 mesh grid the optimum value of w 
is around 1.870. They also point out that when this relaxation parameter is used, the 


convergence is approximately 30 times faster than the usual process (w=1). [Ref. 8] 


10 


II. GOVERNING EQUATIONS 


A. DERIVATION OF THE EQUATIONS 
1. Flow Inside the Circular Tube, Constant Wall Temperature 
Case 


For flow in a circular tube, the viscous energy equation in cylindrical coordinate 


system 1S 
tee 1 ofp H)4 52 pee +202] = 0 (3.1) 
Ox Or |ror\ Or) r°O@\ OG) Ox\ oO 


where di = cdt + aviP and v. is the radial velocity [Ref.1]. Then the equation becomes 


for the constant k 


or ar lk Of at) k @tr . at 
- + es er eee —_— eee Oy = 0 Baz 
ar eS E Or ( =| r* ag" a aa 


The temperature distribution is symmetric (0+/ ag?)=0 and the flow is 


hydrodynamically fully developed, v, = 0. Equation then becomes 


2 
ror\ or) &x k / Ox 


a=k/ pc, the thermal diffusivity of the fluid. So the governing equation for the circular 


duct flow is 


Z 2 
Or? rOr @x? ax 


Now, to non-dimensionalize the problem, assume that 


1] 


{, —T 
6=— r’'= 
fo —ft, 





The governing equation becomes 
2)\)00 076 1 00 #1 066 
l-r =~ ry (3.5) 
Oo Gr’ fF 1On ebecaaey 
2. Flow Inside the Circular Tube, Constant Heat Flux Case 


The governing equation for the circular duct flow is (3.4) 


ot 1a or _ uot 


ar? rer &? ax 
Now, to non-dimensionalize the problem, assume that 


=— = u =— — 
Reyer 


Z 
= , u* =2(1-r* ) 
Qo alk ‘9 Vv 








The governing equation becomes 


| 2\00 076 +1 00 #1 «2876 
l-r — 
Os pr er ncr | Pe? ax* 





(3.6) 


3. Flow Between the Parallel Plates, Constant Wall Temperature 
Case 


The viscous energy equation in cartesian coordinate system for the flow between 


the parallel plates is 


Oi Ol Cl VOX ace Ones 6) GU 
up— +vp— + wp—-| —| k— |+ —] A— 1+ —| k— |} =9 G7) 
Ox Oy OZ (| 00 OG) SOV OV OZ Oz 


where di = cdt + | op and v, is the radial velocity [Ref.1]. Then the equation becomes 
p 
for the constant k 
u ey La 2(2\.2 ee +2(2 = 0 3.8 
Op VS Is edly, | ll Ge Ge Be ey 
The temperature distribution on z direction 1s constant, = = 0 and the flow is laminar and 
z, 
hydrodynamically fully developed, v=0, w=0. Equation then becomes 
2 2 
ede ‘Gr (3.9) 
Ox” Oy k JO 


a=k/ pc, the thermal diffusivity of the fluid. So the governing equation for the parallel 


plates flow is 





2 2 
ce F at _ 4 ot (3.10) 
Ox” Oy” ao 
Now, to non-dimensionalize the problem, assume that 
= pica 8: 2 
Ze Ve ae vie xX un yy as ) 
lye, h V Re Pr 2 
The governing equation becomes 
Z 2 + 
| alate, ee) _u 00 (3.11) 














+ 
Pep, at a 2 Cul 
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4. Flow Between the Parallel Plates, Constant Heat Flux Case 
The governing equation for the parallel plates flow is (3.10) 


6° . aru Ot 


ax? ay a &x 


Now, to non-dimensionalize the problem, assume that 





gq, D,/k h° y° RePr ’ 2 
The governing equation becomes 
| Be-Oeec- ‘ 
_ 5s RICE cE Ee (3.12) 


Pedy ge? gy 2 


5. Effects of Wall Conduction, Flow Inside the Circular Tube, 
Constant Wall Temperature 


The governing equation for the circular duct flow is (3.4) 


6°r lat O7t udt 
—— + —-— + —— = — — 
Gr- rOr @x* ax 


Now, to non-dimensionalize the problem, assume that 


mole py pa (ado u* =2(1-r* ) 
ty —t Yo V Rerr 
The governing equation becomes 
e\d0 09 1 60. 1 86 
is BNR as care Geer as a (3.13) 
Cx Ayn eT or re 
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B. BOUNDARY CONDITIONS 


1. Flow Inside the Circular Tube, Constant Wall Temperature 


Case 


Insulated wall 


d d f f ( Tyan= 19 = constant 


Constant wall temperature 


Insulated wall 





Figure 3. The boundary conditions of circular tube, constant wall temperature flow 


The governing equation is (3.5) 











2)00 076 1 660 #1 866 
]-r = 3 + — ae 3 
Ere GR Gp ea 


where the parameters of the equation are defined above. The boundary conditions for this 








case is: 
oe feleg ie 

at r =0, es 0 (Symmetry boundary condition) 
r 
9 =0 . 

Ate Ay , me x” 20 

= or + 

ore x = 0) 

asx. 3-0, tor, @>1 

asSx° >+0, [t>l, @>0 


sles 


(3.14a) 


(3.14b) 


(3.14c) 


(3.14d) 


2. Flow Inside the Circular Tube, Constant Heat Flux Case 


se///e/h | i dt 


Constant heat flux 





nsulated wall | | | | | 


Figure 4. The boundary conditions of circular tube, constant heat flux flow 


The governing equation is (3.6) 


eee 6°90 106 = #1 «676 
l-r = os en ee oe oe 
Ox Ort or Or Pe at 


where the parameters of the equation are defined above. The boundary conditions for this 

















case is: 
at r =0, = —0 (Symmetry boundary condition) (3.15a) 
r 
00 _itl 
ese ort 2 JO aN (3.15b) 
0g =) for x? <0 
ar* 
asx* 53-0, tor,, G80 (3.15e) 
4 2 
Fae ee ey 
as x* >+0, G@—>-2x" +|—+ - (3.15d) 
48 8 2 
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3. Flow Between the Parallel Plates, Constant Wall Temperature 


Case 


Insulated wall 
a ( / as Tuan = To = constant 


= 3 
uy) . 
3 Constant wall temperature 


Insulated wall 





Figure 5. The boundary conditions of parallel plates, constant wall temperature flow 


The governing equation is (3.11) 


1 66 670 =o u* 86 
Be 5 +4 ; =e 
Epp ax* dy" 














where the parameters of the equation are defined above. The boundary conditions for this 








case 1S: 
at y =0, “ =) | (Symmetry boundary condition) (3.16a) 
6=0 . 
at y"=1, 30 jor x 20 (3.16b) 
dy" = for x <0 
asx*—>-w, [t>r,, @->1 (3.16c) 
asx*>+0, [>%, A 50 (3.16d) 


4. Flow Between the Parallel Plates, Constant Heat Flux Case 





Figure 6. The boundary conditions of parallel plates, constant heat flux flow 


The governing equation is (3.12) 











7 00 80 1 66 1 8 


2 


axt at FOr “Pe eye 
where the parameters of the equation are defined above. The boundary conditions for this 


case 1S: 














at y =0, = =0 (Symmetry boundary condition) (3.17a) 
r 
ole ee: 

at y*=1, Or* - for x* 20 (3.17b) 
00 _, for x* 0 
or* 

asx* 5-0, tr, @>0 (3.17c) 

, Bee | | eee 
aSx* 340, tf, @=-2x 152 es | (3.17d) 
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5. Effects of Wall Conduction, Flow Inside the Circular 


Tube, 
Constant Wall Temperature 


Insulated wall 


Constant wall temperature 


we 
Ge APS IG od Goer roa eS 
ee 2° ees SIS RNS TS pects ie, ne 


Se 





Lae 1= 
Insulated wall 


Figure 7. The boundary conditions of circular tube, constant wall temperature including the circular wall 


The governing equation is (3.5) 

















2)\)00 860 100 #1 3886 
l-r = =p = <a 
Ox ft rf or Pe prt 
The boundary conditions for this case is 
at r =0, = =0 (Symmetry boundary condition) (3.18a) 
Owall a O uid 
atr =1, Coen _ guia (3.18b) 
Or or 
00 
A =0 O < 
atr =l+or, or* for x" s0 (3.18c) 
6-0 jor x 0 
where or is — , ris the thickness of the wall, and ro is the inner radius of the tube 
‘ 
asx* 3-0, to4,, O@->1 (3.18d) 
aS x* 340, [f[>5%, 48-50 (3.18e) 
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IV. NUMERICAL REPRESENTATION 
A. FLOW INSIDE THE CIRCULAR TUBE, CONSTANT WALL 
TEMPERATURE CASE 


1. The Governing Equation 


The governing equation is (3.5) 


Aa oe B a m a Per a2 
Or? r le en ox? 








f1-r?) 28 28 100 1 60 


where 





The finite difference approximations for the first and second derivatives with error of 





order h’ are (figure 2.) 
eet. ; 
08 z (i+1,/) (?-1,7) i o(h: ) (4.1a) 
Ox" (i,7) 2h, 
@..,-6,.. 
08 = yt) ja) +0(h?) (4.1b) 
Or* (i,j) 2h, , 
OPO Fun FF. — U9 oy?) (4.1¢) 
Ox” «,;) hy 
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670 GO ert) 0 eee 

- as UVES sa 1) (i,7) + o(h;) (4.1d) 
Or™ aj h 
(C8) y 





When these approximations are plugged into the governing equation (3.5), the equation 


becomes 

1 Gary t+ Kinny = 2%. | | Mise + Ke j—v 7 2%, | | sey = FD -(-» Gist ~ Mi | _ 
eet cereal | eeamnere mmm mmr (ieee [AS =e) 

Pe?, he h, Yj) 2h, 2h, (4.2) 

After simplifying the expression (4.2), @, ,, becomes 


le vane! 1 vay, 1 
Boras) oy te 14-8) aes - St 1G, | St (4 8 3 
| lage 2h, 2h} “| Pele? 2h, 2h) “SR hy) “PR 2h) (43) 
Se 2 
ie NE 


G, 


2. The Boundary Conditions 


The boundary conditions are: 


a) Atthe Center 





At r=0 = = (Symmetry boundary condition) 
3 
ee: 1 0@ 
When this boundary condition is applied to the governing equation, the element a 
tg 


becomes ~. After applying the L'Hospital's rule, the governing equation for the given 


boundary condition becomes 
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2S (4.4) 
Om or? Pe &* 


Since the temperature distribution at the center on the radial direction is symmetric, 
Gijon = i j-1): 


After applying the finite difference method, the @, ,. at the boundary becomes 


iy) 


fa ee +0 — g 7 
(i+1,/) Pex h? 2h, (Gale) Peth? Zhi, (i,j+1) h* (4 5) 








O44.) ss : 4 
a vy 
Penh, hy 
b) Atthe Wall 
0... =0 
(i,j) i 
atr =, 60 Ute (4.6) 
=i) JOP 
or* 
FOR x’<0 
x'<0, r=, ae = 
Or* 
The governing equation for the given boundary condition becomes 
2 2 
OD eae ae (4.7) 











Ort Pen ax* 
a 00 
Since the wall is insulated, ie 0, O¢ 541) = Oy ;-1) at the boundary. 


After applying the finite difference approximations, the @,, ,, at the boundary becomes 


'J) 


ZS 


] II Z 
: Fist.) P p ie | + OG-4,3) Pe2he i Gorn| 2 


(ca 2 an (4.8) 


as x* >-00, Ths, 7 (4.9) 


ij) 20 (4.10) 


3. The Bulk Mean Temperature 
The mean temperature at a position “1” is (2.8) 


YO 
Lin(x) = Vie Parfait 


YO 
where V is the average velocity and r, is the radius of the tube. When [t,u2ardr iS 
0 


subtracted from both sides and dimensionless paraineter r is substituted inside the 


equation, the equation becomes 
l 
(t,, -t,)= [(¢-t,)u*2r*dr* 
0 


(tn - t0) 


ony (4.11) 





Let’s define the bulk mean temperature, 0, = 


24 


where t. is the inlet temperature (as x—-o) and t, is the constant wall 
temperature.(Fluid exit temperature as x — +00 ) 


Then the bulk mean temperature 0,, becomes 


0 =4f{1-r" p 7 eg ae (4.11b) 


m(x*) 5 
@ is found out from the numerical calculation. After that @,, is calculated by solving the 


equation (4.11) numerically. 


4. The Nusselt Number 


The Nusselt number is (2.3) 


When the heat transfer equations g,=/A(T,-T,,) and qo = , and the 
Flp=R 


dimensionless parameter r are plugged in (2.3), the equation becomes 


(4.12) 





The finite difference approximation for the equation (4.12) using the backward 


difference expression with error of order h’ becomes 


= Se eee en ee 
Nu = = | (ij) : oo aa (4.13) 
m(t) y 


66599 


where “‘j” is the value which corresponds to r =1. 


Zs 


5. Heat Absorption Upto the Entrance 
Assume a differential volume, dV =2zrdr between x=-0o and x=0. The 


energy stored in this control volume is 


dk =c,Atdm (4.14) 
where At is the temperature difference between entrance temperature and the steady state 


temperature after the heating starts, At=¢,,,—-t, and dm is the mass of the control 


volume, dm = pdV .The energy absorbed in the control volume becomes 


07o 
E.psorbed = 20 C, | \(t-te)rdrdx . (4.15) 
-—o 0 
0’7o 
When 27pc, | |t,rdrdx is subtracted from both sides and (t.-t,) is divided to both 
-«o 0 


sides, the equation (4.15) becomes 


r 





0 ‘oO 
_Aabsorbed _ = 4 | [(@-lrdrdx (4.16) 
(t, Ses: -o 0 
where @ = at . As the non-dimensional parameters of x* = ecu and r* =—and 
a Re Pr he 


their first derivatives are substituted into the equation (4.16), the equation becomes 


0] 
cnrared — —27Pe { |@ = l)r*dr* dx" (4.17) 


3 
(t, —t,) PC plo -0 0 
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B. FLOW INSIDE THE CIRCULAR TUBE, CONSTANT HEAT 
FLUX CASE 
1. The Governing Equation 


The governing equation is (3.6) 

















( “ 66 076 106 #1 «06 
l-r = +— + 
ox” art r* Gr* Pes xt 


where 





The finite difference approximations for the first and second derivatives with the error of 


order h? are (4.1a,b,c,d) (figure 2.) 


08 a Feist.) — F545) +o(h?) 
Ox” @,/) 2h 


x 


OF _— O44, js1 is ii, 5-1 me o(h; ) 
or” (i,j) 2h 


Bs 


0°0 Fit.) a Oii-1,3) = 26; 3) +0( h?) 
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2 2 
Ox* Gj) hy 


676 6 eh 6 
é = (i, j+1) ey 1) (7) + o(h;) 
ie 2 

Or (7,7) hy 





a 


When these approximations are plugged into the governing equation (3.6), the equation 


becomes 


l [Fe a "gaa [feet ea ,__ Sen sfc |§. y2 1 — Fy) 
i,j 


° = 
Ped ie hi Yunl dy 2h, | (4.18) 


After simplifying the expression (4.18), 6, ,. becomes 


tj) 





Yen eel 1 1 4] 1 
6.0. yf pao tS - — 14-874 |) a at NG, on) St i'd SO 
‘ la 2h, 2h,) “| eth? 2h, 2h) “ 2hy, | “ 2hy,,} (4-19) 
Gun) = eae 
| 4 4% 
Penh he 
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2. The Boundary Conditions 


The boundary conditions are: 


a) At the Center 


00 


At r=0 rics 0 (Symmetry boundary condition) 








When this boundary condition is applied to the governing equation, the element Me = 
5 


becomes 7 After applying the L'Hospital's rule, the governing equation for the given 


boundary condition becomes 


ed 25 GE (4.20) 
Ox Or* Pe &* 














Since the temperature distribution at the center on the radial direction is symmetric, 











l 


Guia) = Fu,j-n + 
After applying the finite difference method, the 6, ,, at the boundary becomes 
] ] ] ] 4 
(i+1,/) 0 ae +051 J) oo a +O; i+] 7.2 
ae is meee iselenns 27. ans 
he a 5 5 (4.21) 
2 2 aie Ds 
Penh, hy 
b) At the Wall 
Von 20 
ae fale, Piers 
dtr =, — = (4.22) 
Or 
0 Jor 0) 
FOR x" 20, 
M l, a! Pg + 
er IQ j= es 
Or q, a/k a 
BE 
Fe joty ~ Mi,5-1) _ 
2 


2h, 





+ 


so the equation becomes 


The finite difference approximation for this equation is 
When the boundary condition is applied to the governing equation (3.6), the equation 
(4.23) 


becomes 
O 





As the finite difference approximations are substituted in the equation above, the equation 


1 Z l 
+@,.. ceatiel S ee ere 
| wal hy 9 
(4.24) 














becomes 
oo + @ = 
CEI perp? Jo | pe? p2 
O.4,3) ~ 2 D 
——— + ——— 
Per h? hy 
FOR_x <0 
x <0, r=, ae = 0 
or™ 
The governing equation for the given boundary condition becomes 
ao 861 86 
2. 2° 
ar* Pe D axt 
at the boundary. 


= A; 5-4) 
i,j) 


=0, OG jay 


Since the wall is insulated, —— 


l 


] 
ae —— ae 


becomes 
Oe 
1, 
uo \ Per he 
Ds rn Uy: 
Cerse 2 
Been hy 


After applying the finite difference approximations to (4.25), the 6, 
2 
GI) 2 
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(4.25) 


at the boundary 


(4.27) 


d) Asx7> +o 








aS x” >+0, Oi.) > —2x* +{—+ (4.28) 


3. The Bulk Mean Temperature 


The mean temperature at a position “1” is (2.8) 


2 70 
be = |utrdr 


To 
where V is the average velocity and ro is the radius of the tube. When {[t,u2zrdr is 
0 


subtracted from both sides and dimensionless parameter r is substituted inside the 


equation, the equation becomes 
] 
(t,, —t,) = [(¢-t,)u*2r*ar* (4.29) 
0 


Let’s define the bulk mean temperature, 


0, = (. : a, (4.30a) 
a, 


where t, is the inlet temperature (as x — -o ). 


Then the bulk mean temperature @,, becomes 


l 
6, = 4 flier? pear (4.30b) 
0 


S| 


@ is found out from the numerical calculation. After that 6, is calculated by solving the 


equation above numerically. 


4. The Nusselt Number 
The Nusselt number: Nu = = : 


The heat transfer from the wall: q, = h(t oa lye 


The temperature at the position of (x",r"=1): t,. ._.) =t, — 4, a/k 
The average temperature at (x): t, =t, —@,q, d/k 


When these equations are substituted into the Nusselt number equation, Nusselt 


number becomes: 


Nu = (4.31) 





5. Heat Absorption Upto the Entrance 
Assume a differential volume, dV =2zrdr between x=-c and x=Q. The 


energy stored in this control volume is (4.14) 


dk = c,Atdm 
where Aris the temperature difference between entrance temperature and the steady state 


temperature after the heating starts, At=t,,,)—t, and dm is the mass of the control 


volume, dm = pdV .The energy absorbed in the control volume becomes 


a2 


Oo 
Entsorbed = 2%C | [(t—te)rdrdx.. (4.32) 


a 
—-~o 0 


When both sides are divided by q, d/k , the equation (4.32) becomes 





ey ne 
——aesoreed _ = An | [rdrdx C818) 
To PC p Fo /k —o 0 
tems : ‘ 3 x/ Yo + 
where 6 = — . As the non-dimensional parameters of x° = ——— and r* =—and 
Do d/k err ie 


their first derivatives are substituted into the equation (4.33), the equation becomes 


E 01 
——aesorted _ = —AnPe | [O*dr*dx* (4.34) 
Io PC » lo /k -0 0 


C. FLOW BETWEEN THE PARALLEL PLATES, CONSTANT 
WALL TEMPERATURE CASE 

1. The Governing Equation 

The governing equation is defined in equation (3.11) as 


1 076 ae 0°70 u* (06 





Pen, art am “ 2 Ox* 


where 


— 2) 
gator yey pt 2D on ay, 
ty —t, h Re Pr 





The finite difference approximations for the first and second derivatives with error of 


order h’ are (4.1a,b,c,d) (figure 2.) 
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08 isjy ~— ea +0(h?) 


Ox” (i,j) Ap, 
ca: = js — 9G jy + o(h*) 
Cr Ga Ai 


0°0 Gis + a, ~ 2 








= (i,7) a 
iy = — 
(2,7) x 
7, =~ 
O°O _ Oe: say ee Doi, 5-1) 26; 3) 2 
ee ED OD + (2) 
Or* aj) h, 
ig 


When these approximations and wu are plugged into the governing equation (3.11), Fi iy 


becomes 


C] 


ae 








3Yiij BY 4 4 
Chen hee JS oie 7"\ Pez hk? 8h, 8h, RwAN Uh) (4.35) 


2, 8 
ie 


nd 





2. The Boundary Conditions 


The boundary conditions are: 


a) At the Center 





At y=0, j=l, = 0 (Symmetry boundary condition) 
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When this boundary condition is applied to the governing equation, the governing 


equation becomes 











] °0 0°O ~=300 
5 g ae (4.36) 
Pen @x* oy ee 


Since the temperature distribution at the center on the radial direction is symmetric, 
Gijny = Fj 
at the center. 


After applying the finite difference method to (4.36), the @, ,. at the boundary becomes 


l 3 l 3 8 
Gasol eta -2| ae Geso| Feta + Z| i a 








Guay = 5 7 (4.37) 
en -+ ——_ 
Pe*h? h: 
b) At the Wall 
O,  =0 for x 20 
aty =1, ae (4.38) 
oy" 7 for x* <0 
FOR x’<0 
x’ <0, y =l, ue = 0 
oy 
The governing equation for the given boundary condition becomes 
SEWN coe, (4.39) 











2 


Ben ax* Oy~ 
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Since the wall is insulated, ls = 0.0... = 0 eat the boundary. 
aye Gj+1) = %,j-1) ov 


After applying the finite difference approximations, the 0, ,. at the boundary becomes 


i ] 8 
Fist.) ; Pe2h: + As Pe2h |: Bi, 1 E 
ES ne (4.40) 


Oi: 5) = 





e) ASX? -© 


asx*>-0, 0,1 (4.41) 


aS x7 >+0, Cn (4.42) 


3. The Bulk Mean Temperature 


The mean temperature is defined in equation (2.9) as 
l h 
t_ =— jutdy. 
mt Se J ly 


When the dimensionless parameters u’, y are substituted in (2.9), the ty becomes 


iia 
0 
where V is the average velocity and h is the distance from the centerline to the plate. 


When ¢, = |t,uw"dy* is subtracted from both sides, the equation becomes 
2 
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(t -4,)= Ic, —t)utdy* 
0 


Let’s define the bulk mean temperature, 





ae Gate) (4.43a) 


ee) 
where t. is the inlet temperature (as x—-0co) and t is the constant wall 
temperature.(Fluid exit temperature as x > +00) 


‘ y oy Ly . . 


@ is taken from the numerical calculation. After that 6, 1s calculated by solving the 


equation (4.43) numerically. 


4. The Nusselt Number 


The Nusselt number is defined as (2.3) 


When the heat transfer equations q, =i — 1G — peclcs , and the dimensionless 


y=h 


parameter r are plugged in (2.3), the equation becomes 


(4.44) 





a 


The finite difference approximation for the equation (4.44) using the backward difference 


expression with error of order h” becomes 


4 ( 30; , ~40;; 1) +0; ;- 
a = mee | (4.45) 
a 


my 


5. Heat Absorption Upto the Entrance 
Assume a differential volume, dV = bdxdy between x=-0co and x=0. The 
energy stored 1n this control volume is (4.14) 
dE =c,Atdm 
where Af is the temperature difference between entrance temperature and the steady state 


temperature after the heating starts, At=%,,,—¢, and dm is the mass of the control 


volume, dm = pdV .The energy absorbed in the control volume becomes 
0 +h 


| Deen ee = bpc | {(¢ —t,)dydx a (4.46) 


—-xo-h 


0h | 
When 2dpc, | |t,dvdx is subtracted from both sides and (t, -t,) is divided to both 


—o 0 
sides, the equation (4.46) becomes 


0h 
Len aon -9 | | (0 = 1) dydx (4.47) 


bpc,(t, =) ~o 0 
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where 0@= Geom . As the non-dimensional parameters of x* = elses ee 
(Z, ins to) Re Pr 


r~ =—and their first derivatives are substituted into the equation (4.47), the equation 


lo 
becomes 
absorbed _ _ —4Pe { (@ —l)dy*dx* (4.48) 
h*bpe(t, —t,) —co 0 | 


D. FLOW BETWEEN THE PARALLEL PLATES, CONSTANT 
HEAT FLUX CASE 
1. The Governing Equation 
The governing equation is (3.12) 
1 06 ,0°6 u* 06 
at eS 
Poa at a 2 Ox* 


where 


tt 4p ee, 


te ., , Xx ,  Dp=4h 
q. D, /k h Re Pr 





The finite difference approximations for the first and second derivatives are (4.1a,b,c,d) 


06 Fist Fay +o0(h2) 
Ons) 2h, 
a _ Fein ~ 85-0 oh?) 
Or” (i,j) 2h, 
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O70 Minty t ayy ~ 29%, +0(h?) 








5 Z 
Ox” «,7) Nx 
2 = 
OO OG jay t 9 jy — 29,59 2 
= a + (hh) 
ar”? h; ; 
(J) y 


When these approximations and u’ are plugged into the governing equation, the bulk 


temperature at the position (i .j) is written as 





ee: 1 23ye 4 4 
Gosol BEE 2p? a a 8h. +0641, Petia rane +95 ja) a) +O, 5-1) 7) 
Pal Sed ke Sees ak ENE EE (4.49) 
Qj) 2 g 





+ — 
ee = 
Pe2h? 


2. The Boundary Conditions 


The boundary conditions are: 


a) At the Center 





At y=0, j=l, =0 (Symmetry boundary condition) 


When this boundary condition is applied to the governing equation (3.12), the governing 
equation becomes 
1 6@ hs 6°60 3 60 


2 ra ao oe a (4.50) 
Pen ax* Oy 4o%x 











Since the temperature distribution at the center on the radial direction is symmetric, 


A. ja) = 9G. ;-1) at the center. 
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After applying the finite difference method, the 0, ,. at the boundary becomes 


G.. a +O), Se HOG 4 
vas Pe2 2 2] ( ol sei 2] ue {5 4 








Fea) = a: 
——— — +--—. 
Penh? h* 
b) At the Wall 
0 OF ae 0 
a 00 J 
at y =1, —=) ] 
4 US 0) 
FOR_x <0 
me) ll, a =0 
Oy 


The governing equation for the given boundary condition becomes 


zl ACR, (4.52) 
Pen ax* dy* 





Since the wall is insulated, ee =o eo — 0 eat the DOUNGaLy, 
Gy* (7+) @,J-1) 


After applying the finite difference approximations, the @,, ,. at the boundary becomes 


] J 8 
Oe: a. +6, st | oreeemen oa + 0, are 


tee 
Pe*h? hy 


Ou) = (4.53) 


4] 























FORM 22 0 
Sa = 0 ‘ y =] : OF = —— 
oA 
The governing equation for the given boundary condition becomes 
LO eG (4.54) 
ee Ou Oy* 
Cs oe h 
0g = _ Giy+l) CUE), = i : Oi; ja) —_— 0; j-) itt de at the boundary. 
Oy” + Piss + 2 
After applying the finite difference approximations, the 6, ,. at the boundary becomes 
J ] 8 2 
sul esa) Sool pas) Oo) a 
= an (4.55) 
By 9 Be 2) 
LAO en Oe 
Cc) Asx 7 -© 
as x* 3-0, O., 29 (4.56) 
d) Asx > +o 
39 3 ae 2 
as x* O,, 5 =—2x7 +| —— + — = Ae 
a GD) 1120 16 6 a 


3. The Bulk Mean Temperature 


The mean temperature is (2.9) 


] ? 
= — |utdy. 
mn ee 


es 
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When the dimensionless parameters u', y are substituted in (2.9), the tm, becomes 


l 
t, = lu Tdy” 
0 
where V is the average velocity and h is the distance from the centerline to the plate. 


When ¢t, = |t,u*dy” is subtracted from both sides, the equation becomes 
0 


l 
(t,, -t,)= IG —t, judy” 
0 
Let’s define the bulk mean temperature, 


= Sete), (4.58a) 
och 


m 


where t, is the inlet temperature (as x + -oo). Then the bulk mean temperature 6,, 


becomes 


1 
0 = ; J f z y” bay’ (4.58b) 


@ is taken from the numerical calculation. After that 0,, is calculated by solving the 


equation (4.58) numerically. 


4. The Nusselt Number 


The Nusselt number: Nu = D, 





The heat transfer from the wall: g, = Ce eo 
The temperature at the position of (x",r'=1): ¢.. .) =t.- 6q, D,, /k 
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The average temperature at (x): t, =, —9,,q, D,/k 
When these equations are substituted into the Nusselt number equation, Nusselt number 


is Calculated as 


(4.59) 





5. Heat Absorption Upto the Entrance 
Assume a differential volume, dV = bdxdy between x=-00 and x=0. The 
energy stored in this control volume is (4.14) 
dk =c,Atdm 
where Afis the temperature difference between entrance temperature and the steady state 


temperature after the heating starts, At=%, ,,)—f, and dm is the mass of the control 


volume, dm = pdV .The energy absorbed in the control volume becomes 


0 +h | 
Le heofhed = bpc, | \¢ —t, )bdydx : (4.60) 
h 


—-mO — 


When both sides are divided by q, D,,/k , the equation (4.60) becomes 





ST peorbed a 
— absorbed _=_8 | [Odydx (4.61) 
FPLC, h/k -o 0 
os ‘ : 
where 0 = —4 . As the non-dimensional parameters of x” = Bes and r* =——and 
Do D,,/ lo 


their first derivatives are substituted into the equation (4.61), the equation becomes 


IO Pe: d Pa es ea 
—sosorore = —16P Ody~ dx 4.62 
7 bit /a “t ! ly (4.62) 
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E. THE EFFECTS OF WALL CONDUCTION, FLOW INSIDE THE 


CIRCULAR TUBE, CONSTANT WALL TEMPERATURE CASE 


1. The Governing Equation 


The governing equation is (3.5) 














fir?) 28 £8 100 1 86 


where 








The finite difference approximations for the first and second derivatives with error of 


order h? are (4.1a,b,c,d) (figure 2.) 





a0 OO. 
ea 
Ox (i,J) PNG 
OO 6; +1 —6;; —l 
- = (?,j/+1) Cy, ) + o(h:) 
Or Cy) on, 
6°6 Caan, 7 Oii-4. 4) Ee 20); 5, 2 
= = i aes + o(h; ) 
Ox” i,j) x 
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6°64 0, yee any = 2 
_ (i, j+1) (i,j-1) (i,j) + o(h>) 


2 2 
Ch ac h, 





When these approximations ar: plugged into the governing equation (3.5), the equation 


becomes 


I GKissj) + inj — 24,9) | z Sue + OG, j-1) ~ 28 4i,5) | __! en ~ Gu, j-1) } ( a [fo = Hi-1,3) 


3 » Tod] = 
Per h, hy Mi.j) 2hy i a | (4.63) 


After simplifying the expression (4.63), the equation becomes 


ree ee leon l l 1 ! 
Boas) may te 1 ery) os Ot 1 | StF 8 -)| GS - 
vol Fa 2h, 2h) “| Peek? 2h. 2h) “A thy, | “A? hy) (4-04) 
oe. a aa ears: imei’ =< Sane a ee 
——_—_—_- + — 
Pesh: he 


O14.) = 


2. The Boundary Conditions 


The boundary conditions are: 


a) At the Center 


At r=0, =0 (Symmetry boundary condition) 





00 


+ 





When this boundary condition is applied to the governing equation, the element a 
ron 


becomes 7 After applying the L'Hospital's rule, the governing equation for the given 


boundary condition becomes 


08 _,00 | 1 oé (4.65) 


ex? Be Pe’ Ax’ 











Since the temperature distribution at the center on the radial direction is symmetric, 
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0 0 


Gyo)  %Gj-)° 


After applying the finite difference method, the 6,, ,. at the boundary becomes 


l l l : 
6. ae +6. SaaS oe spare 2. 
wl Sep +] Gi ol sa Zin | GJ fs | 
Ou 3) = i een co =| | ae Se — 
Pephy hy 
b) At the Wall 
Owall = 9 uid 
at r mall: 061,01 = CF uid 
art or* 








86,0; OF uid 
K solid = =k fluid are 


Using the 2"° order forward and backward difference method, the equality 


becomes 


30, 5 — 40; 4) + &; ;- —36,, , + 4O,; a) — Fi; | 
if G7) = 1) (i,j |. | (i,j) aa ee) (4.67) 
uf ys 


If hys=hye, the bulk temperature at the fluid-wall boundary becomes 


k, k, 
4 Gi, jay + F~%i,5- || %es42) + F,3-) 


k 
c + Z| 


Oi) = (4.68) 
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Cc) On the Outer Edge of the Wall 
FOR x <0 


For x <0 and r* =1+6y,_ the outer surface of the solid wall is insulated where 


_ solid _wall _ thickness 
tube _inner_ radius — 
The governing equation for the given boundary condition becomes 


2 2 
se (4.69) 
r ae ee” 








Since the wall is insulated, = =0,0@4 =@ _1 at the boundary. 
r 


2 eg F1) (J 


After applying the finite difference approximations, the 0 


(i,j 


] ] Z 
am Pez he | a am Pe2h + OG 5-1 2. 


, at the boundary becomes 


O04) = 5 5 (4.70) 
ee + —— 
Pezh, hy 
FOR x >1 
Oi: ;) — 0) (4.71) 
d) ASX? - a 
as x* >-0, Ga, 1 (4.72) 
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asx*>+0, 6,0 (4.73) 


3. The Bulk Mean Temperature 


The mean temperature at a position “1” is (2.8) 


YO 
L(x) = Ve Jecsosteanr@r 


TO 
where V is the average velocity and r, is the radius of the tube. When [¢,w2zrdr is 
0 


subtracted from both sides and dimensionless parameter r is substituted inside the 


equation, the equation becomes 


] 
(t,, —t,)= {(t-t,)u* 2r*dr* 
0 


Ga) 


Let’s define the bulk mean temperature, 0, = oy where t, is the inlet temperature 
e “o 


(as x —> —co) and t, is the constant wall temperature.(Fluid exit temperature as x — +00) 


Then the bulk mean temperature 0,, becomes 
ee =4f-r? ‘dr’ (4.74) 
0 


@ is found out from the numerical calculation. After that ,, is calculated by solving the 


equation (4.74) numerically. 
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4. The Nusselt Number 


The Nusselt number is defined as (2.3) 


_hd 


Nu 
k 


When the heat transfer equations g,=h(t,—-t,,) and qo = , and the 
F lpsR 


dimensionless parameter r' are plugged in (2.3), the equation becomes 


(4.75) 





The finite difference approximation for the equation (4.75) using the backward 


difference expression with error of order h? becomes 


2 { -36, , +46, 4, -9%, j> 
Nie a F . Cap 1) (i,j ; (4.76) 
m(i) uy 


where “j” is the value which corresponds to r =1. 


5. Heat Absorption Upto the Entrance 
Assume a differential volume, dV =2zrdr between x= -0o and x=0Q. The 


energy stored in this control volume is (4.14) 


dk = c,Atdm 
where Aris the temperature difference between entrance temperature and the steady state 


temperature after the heating starts, At=1,,)—t, and dm is the mass of the control 


volume, dm = pdV .The energy absorbed in the control volume becomes 
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079 
Expsorbed = 270, | [(t—te)rdrdx . (4.77) 


a 
-x 0 


0 7o 
When 27pc, { |t,rdrdx is subtracted from both sides and (t--t.) is divided to both 
—-~o 0 


sides, the equation (4.77) becomes 


Oo 





“absorbed = 24 { [(@-Drdrdx (4.78) 
(t, at o) pe p —-o 0 
t—t, : . sages 2) 00: 2 
where 0 = . As the non-dimensional parameters of x” =——— and = — and 
Ln =, Re Pr i 
their first derivatives are substituted into the equation (4.16), the equation becomes 
Ho hsorhed — —J7P ar + + dx* 
ae ae ee |@-Dr*dr (4.79) 


« =i one —0 0 


3) 





V. NUMERICAL RESULTS AND DISCUSSION 


A. BULK MEAN TEMPERATURE VS. X 


The bulk mean temperature graphs of the fully developed laminar flow inside the 
circular tube and parallel plates for constant wall temperature and heat flux cases were 
plotted in figures 8-11. The Peclet numbers 0.5, 1, 2, 3, 5, 10, 50, and 100 were used for 
the plots to see how the Peclet number, Pe=RePr, affects the axial conduction in the flow. 
As it is seen from the figures, when the Peclet number gets smaller, the axial conduction 
effects of the fluid becomes very important and the temperature of the fluid upto the 
entrance increases. The axial conduction of the fluid has to be considered for these cases 
and must not be neglected for the Peclet numbers less than about 10. When the Peclet 
number gets larger, the effects of the axial conduction decrease. As it is seen from the 
figures 8-11, for the Peclet numbers 50 and 100, the bulk mean temperature of the fluid at 
the position x*=0 is almost equal to the entrance temperature, which means that there is 
almost no heat absorption or temperature increase upto the entrance. The axial conduction 
of the fluids can be neglected for these high Peclet numbers. 


1. Flow Inside the Circular Tube, Constant Wall Temperature 


Case 


The graphs of the bulk mean temperature vs. x for the circular tube, constant wall 


temperature flow was plotted for various Peclet numbers in figure 8.a and 8.b where theta 


bulk meanis 6, = 


The validity of the figure 8.b has been checked by comparison with Ref.5. 


Theta Bulk Mean 


CIRCULAR TUBE CONSTANT WALL TEMP. FLOW 





» THETA BULK MEAN 
TN Sasa each ree so cams aces val 1-155 


and the nondimensional axial distance x” is x* = 





2a Ga. 
Re Pr 


Figure 8.a The axial bulk mean temperature distribution for circular tube constant wall temperature flow 


Theta Bulk Mean 


CIRCULAR TUBE CONSTANT WALL TEMP. FLOW, THETA BULK MEAN 


1 = I RPO OO et 
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OE) a i a Be ee eo 


Oh ee Se a ee ee ine Sa ae 


OFT ee ee ee ee ee hh = 


Figure 8.b The enlarged view of figure 8.a 
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Pe =0.5 
Pe=1 
Pe=2 
Pe=3 
Pe=5 
Pe =10 
Pe =50 
Pe =100 











2. Flow Inside the Circular Tube, Constant Heat Flux Case 


The graphs of the bulk mean temperature vs. x for the circular tube, constant heat 


flux flow was plotted for various Peclet numbers in figure 9.a and 9.b where theta bulk 


mean is @,, = (0. ='m) and the nondimensional axial distance x” is x* = cue , 
q, a/k Reva 


CIRCULAR TUBE CONSTANT HEAT FLUX FLOW, THETA BULK MEAN 


Theta Bulk Mean 





Figure 9.a The axial bulk mean temperature distribution for circular tube constant heat flux flow 


CIRCULAR TUBE CONSTANT HEAT FLUX FLOW, THETA BULK MEAN 


Theta Bulk Mean 
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Figure 9.b The enlarged view of figure 9.a 
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3. Flow Between the Parallel Plates, Constant Wall Temperature 


Case 


The graphs of the bulk mean temperature vs. x” for the parallel plates, constant wall 


temperature flow was plotted for various Peclet numbers in figure 10.a and 10.b where 


theta bulk meanis g = ea and the nondimensional axial distance x is x+ = 2x1 Dy | 
"  (t, -t, Re Pr 


Theta Bulk Mean 





0 
-20 -15 -10 5 0 = 


Figure 10.a The axial bulk mean temperature distribution for parallel plates constant wall temperature flow 


PARALLEL PLATES CONSTANT WALL TEMP. FLOW, THETA BULK MEAN 
q , 


= re ee ee eae 





Theta Bulk Mean 


Figure 10.b The enlarged view of figure 10.a 
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4. Flow Between the Parallel Plates, Constant Heat Flux Case 


The graphs of the bulk mean temperature vs. x for the parallel plates, constant heat 


flux flow was plotted for various Peclet numbers in figure 11.a and 11.b where theta bulk 


oasioae | 
mean is @,, = (e=tn) and the nondimensional axial distance x” is x* = ea 
Go Dy, /k Re Pr 


PARALLEL PLATES CONSTANT HEAT FLUX FLOW, THETA BULK MEAN 


Theta Bulk Mean 


Theta Bulk Mean 





Figure 11.b The enlarged view of figure 11.a 


B. NUSSELT NUMBER VS. X 


The nusselt number graph of the fully developed laminar flow inside the circular 
tube and parallel plates for constant wall temperature and heat flux cases were plotted in 
figures 12-15. The Peclet numbers 0.5, 1, 2, 3, 5, 10, 50, and 100 were used for the plots 
to see how the Peclet number, Pe=RePr, affects the nusselt number distribution in the 
flow. 


1. Flow Inside the Circular Tube, Constant Wall Temperature 


Case 
The graph of the nusselt number vs. x for the circular tube, constant wall 
temperature flow was plotted for various Peclet numbers in figure 12 where x” is 


5 Opeliie 
i 
Re Pr 





CIRCULAR TUBE CONSTANT WALL TEMP. FLOW, NUSSELT NUMBER 


Nusselt Number 





Figure 12. The nusselt number distribution for circular tube constant wall temperature flow 
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The nusselt numbers approach to the values in those of table 1 when x gets larger as seen 


in figure 12. 
Pe 0.5 I Z 3 5 10 50 100 
Nu. no | 4.0971 4.0276 3.9224 3.8506 3.767 3.6948 3.6586 3.6572 


Table 1. The nusselt number for various Peclet numbers for circular tube constant wall temperature case as 
x goes to infinity 


2. Flow Inside the Circular Tube, Constant Heat Flux Case 


The graph of the nusselt number vs. x for the circular tube, constant heat flux flow 


2x/d 
Re, Pr 





was plotted for various Peclet numbers in figure 13 where x° is x" = 


CIRCULAR TUBE CONSTANT WALL TEMP. FLOW, NUSSELT NUMBER 


Nusselt Number 





Figure 13. The nusselt number distribution for circular tube constant heat flux flow 


The nusselt number is independent from the Peclet number for this case, and approaches 


to 4.364 where x goes to infinity. 
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3. Flow Between the Parallel Plates, Constant Wall Temperature 


Case 


The graph of the nusselt number vs. x for the parallel plates, constant wall 


temperature flow was plotted for various Peclet numbers in figure 14 where x” is 


/D 
et cai led 5 is 4h. 
Re Pr 


PARALLEL PLATES CONSTANT WALL TEMP. FLOW, NUSSELT NUMBER 


Nusselt Number 





Figure 14. The nusselt number distribution for parallel plates constant wall temperature flow 


The nusselt numbers approach to the values in those of table 2 when x gets larger as seen 
in figure 12. 

Pe 0.5 I 4 3 o 10 50 100 
Nu.no | 8.0588 8.0059 7.9168 7.8464 7.7471 7.6303 7.5457 7.542 


Table 2. The nusselt number for various Peclet numbers for parallel plates constant wall temperature case 
as x goes to infinity 
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es 
x Flow Between the Parallel Plates, Constant Heat Flux Case 


The graph of the nusselt number vs. x for the parallel plates, constant heat flux 


2x/D,, 
Ke EF 





flow was plotted for various Peclet numbers in figure 15 where x” is x” = , D, 1S 


4h. 


PARALLEL PLATES CONSTANT HEAT FLUX FLOW, NUSSELT NUMBER 
3 


Nusselt Number 





Figure 15. The nusselt number distribution for parallel plates constant heat flux flow 


The nusselt number is independent from the Peclet number for this case, and 


approaches to 8.235 where x goes to infinity. 
C. HEAT ABSORPTION UPTO THE ENTRANCE 


The energy change from the entrance to x=0 where the heating starts has been 
investigated and plotted for fully developed laminar flow inside the circular tube and parallel 
plates for constant wall temperature and heat flux cases in figures 16-19. The Peclet numbers 
0.5, 1, 2, 3, 5, 10, 50, and 100 were used for the plots to see how the Peclet number, Pe=RePr, 


affects the heat absorption of the fluid. As it is seen from the graphs, the heat absorption 
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logarithmically increases as the Peclet number decreases. Therefore axial conduction effects or 
heat absorption upto the entrance where the heating starts has to be considered for low Peclet 


numbers. 
1. Flow Inside the Circular Tube, Constant Wall Temperature 
Case 


The graph of the heat absorption, E’ vs. Peclet number for the circular tube, 


constant wall temperature flow was plotted in figure 16 where the dimensionless 


absorbed energy E is p* _ __ Eabsorbed 
3 
(t, =i Wear 


, ENERGY ABSORPTION, CIRCULAR TUBE, CONSTANT WALL TEMP. 


0 upto the Entrance, E* 


Energy Absorbed from x 





Peclet Number 


Figure 16. The energy absoption upto the entrance for circular tube constant wall temperature flow 


2. Flow Inside the Circular Tube, Constant Heat Flux Case 


The graph of the heat absorption, E’ vs. Peclet number for the circular tube, 


constant heat flux flow was plotted in figure 17 where the dimensionless absorbed energy 


eee E 
E is E _. ___ absorbed __ 
gn hey ta [ke 
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ENERGY ABSORPTION, CIRCULAR TUBE, CONSTANT HEAT FLUX 


0 upto the Entrance, E° 


Energy Absorbed from x 





Peclet Number 


Figure 17. The energy absoption upto the entrance for circular tube constant heat flux flow 


3. Flow Between the Parallel Plates, Constant Wall Temperature 


Case 


The graph of the heat absorption, E* vs. Peclet number for the parallel plates, constant 


wall temperature flow was plotted in figure 18 where the dimensionless absorbed energy E” is 


E’ =— E absorbed ; 
h’bpe,(t, -te) 


: ENERGY ABSORPTION, PARALLEL PLATES, CONSTANT WALL TEMP. 
10 


0 upto the Entrance, E° 


-”~ 
= 


Energy Absorbed from x 





Peclet Number 


Figure 18. The energy absoption upto the entrance for parallel plates constant wall temperature flow 


4. Flow Between the Parallel Plates, Constant Heat Flux Case 


The graph of the heat absorption, E* vs. Peclet number for the parallel plates, constant 


heat flux flow was plotted in figure 19 where the dimensionless absorbed energy E’ is 


E a E absorbed 


bh’ /a 


; ENERGY ABSORPTION, PARALLEL PLATES, CONSTANT HEAT FLUX 
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0 upto the Entrance, E* 


Energy Absorbed from x 





Peciet Number 


Figure 19. The energy absorption upto the entrance for parallel plates constant heat flux flow 


D. AXIAL CONDUCTION EFFECTS INCLUDING THE WALL 


CONDUCTION 


The axial conduction effects for circular tube constant wall temperature case has 
been investigated including the wall conduction effects. The Peclet numbers 1, 3, 5, 
outer/inner radius ratios of 1.2, 2 and fluid/wall thermal conductivity ratios of 1 and 1/10 
were used to see how these properties affect the temperature, nusselt number and heat 


absorption. The theta bulk mean temperature, nusselt number variations were plotted vs. 
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x” and the heat absorption were quantified as a function of Peclet number, thermal 
conductivity and radius ratios. The results are seen on the following pages. 


1. Bulk Mean Temperature Variations 


The _%:_ vs. x” graphs are seen in figures 20.a and 20.b for different Peclet 


numbers, thermal conductivity and wall thickness. Since the bulk mean temperature 
values are very small about x*=0, the normalized bulk mean temperature on the centerline 


is a good way to see the change at the bulk temperature relative to the centerline 


(tte) aad x ig xt = 2 
(= =2,) Re Pr 





temperature at x =0, where the bulk temperature is 0 = 


The bulk mean temperature vs. x graphs for the Peclet numbers 1,3, and 5, and the 


thermal conductivity ratios of 1, 1/10 and the radius ratios of 1.2 and 2 are seen in figure 


Game 


21 where the bulk mean temperature is 0, = am. 
aU 


THE NORMALIZED BULK MEAN TEMPERATURE 


Pe=1k=1 h=1.2 
Pe=1 k=1/10 h=1.2 
Pe=3 k=1 h=1.2 
Pe=3 k=1/10 h=1.2 
Pe=5k=1 h=1.2 
Pe=5 k=1/10 h=1.2 
Pe=1 without wall 
Pe=3 without wall 
Pe=5 without wall 





Figure 20.a The normalized bulk mean temperature with and without wall 
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THE NORMALIZED BULK MEAN TEMPERATURE 


Pe=ik=1 h=2 
Pe=3 k=1 h=2 
Pe=5 k=1 h=2 
Pe=1 without wall 
Pe=3 without wall 
Pe=5 without wall 





Figure 20.b The normalized energy absorption upto the entrance for parallel plates constant heat flux flow 


The theta bulk mean values with the wall conduction case is closer to the steady 
state values than the no wall conduction case as it is seen in figure 21. The bulk mean 
temperature go to zero as x goes to plus infinity. This means that the temperature of the 
fluid including the wall conduction increases faster than the temperature of the fluid 
neglecting the wall conduction. The main reason of the relative increase at the 
temperature including the wall conduction case is the heat that is transferred by the axial 
conduction of the wall opposite to the flow direction. The heat moves to the opposite of 
the flow direction and is absorbed by the fluid. So the fluid temperature increases much 
more than no wall case. 


THE BULK MEAN TEMPERATURE VS. X 


Pe=1 k=1 h=1.2 
Pe=3 k=1 h=1.2 
Pe=5 k=1 h=1.2 
Pe=1 k=1/10 h=1.2 
Pe=3 k=1/10 h=1.2 
Pe=5 k=110 h=1.2 
Pe=1 k=1 h=2 
Pe=3 k=1 h=2 
Pe=5 k=1 h=2 
Pe=1 without wall 


Pe=3 without wall 
Pe=5 without wall 





Figure 21. The bulk mean temperature with and without wall conduction 
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2. Nusselt Number vs. X 
The Nusselt numbers vs. x+ graphs are seen in figure 22 for different Peclet numbers, 


2xia 
Re Pr 


. The Nusselt number 





thermal conductivity and wall thickness where x" is x* = 


distribution is very different from the no wall conduction case as it is seen in the figure. 
As the wall gets thinner and the thermal conductivity increases, the values of the Nusselt 
number gets closer to the case of no wall case. h=1/10 and k=1.2 1s close to the no wall 


conduction case, so the values of the Nusselt number for these curves are close to the 


values of no wall case. 


THE NUSSELT NUMBER VS. X 
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Figure 22. The nusselt number distribution including the wall conduction 
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3. Heat Absorption vs. X 


The graphs of the heat absorption, E* vs. Peclet numbers for the conductivity ratios of 


1,1/10 and radius ratios of 1.2 and 2 were plotted in figure 23 where the dimensionless absorbed 


energy E” is f° = ——Pabsorted __. As it is seen from the graph, As the wall's conductivity or 
(eteece 


thickness increases, the heat absorption before x=0 significantly increases. The wall 
conduction case has to be added to the problem for the microheat exchanger problems. 
Because the relative thickness of the wall is very large as it is compared to the inner 
diameter of the tube. The values of h=1.2 or h=2 are realistic for the microheat machines, 
because the diameter of the tubes are very small for these heat exchanger. The heat 
absorption for kf/ks=1/10 and h=1.2 case is more than ten times of the case of no wall. If 


the wall thickness would be h=2, then the heat absorption would be much greater. 


ENERGY ABSORPTION VS. X 


/10 h=1.2 
oO wall case 


Peclet Number 





Figure 23. The energy absorption upto the entrance with and without the axial wal! conduction 
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VI. CONCLUSIONS 


Axial conduction effects in duct flows are important and have to be considered for 
cases of Pe<l10. The theta bulk mean temperature and Nusselt number variations vs. 
downstream distance were plotted as seen in the results. It 1s obviously seen the increase 
in the temperature against the flow direction from x=0 as the Peclet number decreases. 
Heat absorption by the fluid in the insulated region was also quantified by using the non- 
dimensional parameter E’. The heat absorption increases logarithmically as the Peclet 
number decreases, as seen in the heat absorption graphs. When the wall conduction 
effects are included in the problem, it was calculated that most of the heat absorption by 
the fluid is because of the wall conductivity as the wall conductivity gets larger. The same 
effect is also seen with increasing wall thickness. The heat absorption by the fluid gets 
larger when the wall gets thicker. It is very important to include the wall conduction in 
the problem as the wall gets thicker and the thermal conductivity of the wall gets larger. 
With large wall conductivity and wall thickness, the fluid temperature increases 
significantly before the heating starts. The current interests in axial conduction effects are 
because of the advent of micro heat exchanger technology with microtubes where the 


Reynolds number is very small. 
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APPENDIX —- A SAMPLE PROGRAM 


The Program for Circular Tube, Constant Wall Temperature Case 


This program solves the finite difference equations for the circular tube, constant 
wall temperature case by using the Gauss Seidel] Iteration Method with the overrelaxation 
parameter w=1.8. The program solves the equations doing the following steps : 

1. The matrix is formed by the program, the first row represents the region which 
touches the tube. The last row is the centerline. The first column is the condition, where x 
goes to minus infinity. The last column represents the boundary condition as x goes to 
plus infinity. 

2. The boundary conditions are applied to the problem. The centerline boundary 
condition, the boundary conditions at r=r,, the boundary conditions as x goes to minus 
and plus infinity are applied to the last row, first row, first column and the last column 
respectively. 

3. The iterations are made calculating the values at all the grid points. As it is seen 
in the program, the odd rows are calculated first, the even rows later. The convergence of 
the problem is better in this case. 

4. The difference between the error and the boundary conditions for the plus and 
minus infinity is checked. If the boundary conditions and the difference between the 
iterations are reasonable, the program is stopped and the data 1s processed. The max. 
difference between the iterations is less than 1E-6. The max difference between the 
values of the horizontal grid points which represents minus and plus infinity is less than 


1E-6. 
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WF off oP oP ol? 


oe 


oe 


oe 


oe 


oe 


This problem solves the problem for circular tube 


constant wall temperature case 


Programmed by Ibrahim Girgin 


format long 


clear 
Start=1; 
iteration _number=300000; 


ae ole 


oe 


Nx=2501 
Ny=101 
KicsEeec= 501 


% For Peclet number is 1 


Pe=1 


6 For the relaxation parameter 1. 


w=1.8 
save _iteration=100 


8 


Nx is the column number of the matrix 
Ny is the row number of the matrix 
x zero is the column number which represents x=0 


hy 1s the vertical distance between the grids 
hx is the horizontal distance between the grids 


hy=1/ (Ny-1) ; 
hx=hy; 


teta is the matrix which holds the grid points 


initialization of the teta 


teta=ones (Ny, Nx) ; 


Apply the boundary condition for r=Ro, 


p49 


teta(:,x_zero:Nx)=zeros (Ny, length(x_zero:Nx) ); 


6 €.9. a =Ny SS) 
% the initial 5 Ea 
% value of ge os gee 
% theta te eb 
% ike See 1 
ro J ee ee 
% i=l] 


ae ol? 


o\? 


PRR PP 


2 


ee eee 


PRP RPP PB 


The values for this Eheta matrix is 


ee 
ith mole} ro) (=) 
200000 


=<. 2eu° 


oO 2 © © o@ © 


OO © -Oo7o © 


Oo © Oo 2 2 oo 


Oo DO Ooo 


OS Ov O- O20 > 


a 
It 


Nx 


Nx=14 Ny=6 x_zero=8 hx=0..2 hy=0.2 
Pe,w,file name,save iteration is what is to be given 


of oP 


+ The variables that are being used inside the loop 


varl=Pe”2*hx*2; 
var2=2/hy°2; 
var3=1/2/hx; 
var4=1/hy~2; 
varS=1/varl-var3; 
var6=1/varl+var3; 
eee = 1) 27 by; 
var8=2/varl1l+4*varé4; 
var9=2/varl+var2; 
varl0=2*hx; 


o\? 


The iteration starts here. For every value of kkk, one 
% iteration is made 


for kkk=start:iteration_number 


perpintax err=Leta; 


for 12=2:Nx-1 6 do the iteration for all the columns except the 
oer S 
¢ and the last columns 


Mor ie Ls 2 yy % do the iteration for all the odd rows 
y=(g- 1) by; o calculate the y, a= ee=1 
196 3) SS % do the iteration for radius=0 (for the center) 


teta(Ny,1)=(teta (Ny,i+1) * (var5) +teta (Ny,1-1) * (var6) + 
4*teta (Ny-1,i)*var4) /var8; 


elseif j==Ny % do the iteration for radius=R 
ie ere 6 DR a<Ke Zeno & for Eadiucer, st sts Iniculatcd 
teta(1l,i)= ((teta(1l,1+1)+teta(1,i- 


1))/(varl) +teta(2,i)*var2) /var9; 


end 


9 


else % do the iteration for all the odd rows except r=0 
and r=R 


teta (Ny-j+1,i) =teta(Ny-j+1,i) +w* ((teta(Ny-j+1,i1+1)*(var5 + 
ie ay i g-\ a 116) emer 
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+teta (Ny-j+1l,i-1)*(var6 - y°2/varl10)+teta(Ny-j,i)*(var4 + 
var7/y)+... 


teta (Ny-j+2,1)*(var4 - var7/y))/var9-teta(Ny-j+1,i)); 


end 
eng 
EOr =e Ny) % do the iteration for all the even rows 
y=(j-1) *hy; 


teta (Ny-j+1,1) =teta (Ny-j+1,1) +w* ((teta (Ny-j+1,1+1)*(var5 + 
y 2/vari0) 2... 


+teta (Ny-j+1,i-1)*(var6 - y°2/varl0)+teta(Ny-j,i)*(var4 + var7/y) + 


teta (Ny-j+2,i)*(var4 - var7/y))/var9-teta(Ny-j+1,1i)); 


end 


o\? 


end End of the iteration 


oo 


Save the variables at every "Save_iteration” 
e.g. 1£ the save ateratien equals to 20,7 then save 
the variables at when kkk=20,40,60,80,... So on. 


o\? 


oo 


if fix(kkk/save_iteration) ==kkk/save_iteration; 
max error= (max (max (abs (teta-for_max_err)))) 


eleat JEGreMasGic 5 


save forpel teta Nx Ny x_zero w Pe hx hy max_error kkk 
save iteration start y 


end 


end 
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